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^5 Abstract 

We calculate the dominant loop corrections from both standard and supersymmet- 
Qh. ric QCD to the effective coupling of neutralinos to nucleons. The potentially largest 

corrections come from gluino-squark loop contributions to the Higgs boson cou- 
plings to quarks; these corrections also affect the leading spin-independent squark 
exchange contribution. Ordinary QCD corrections to the effective coupling of CP- 
even Higgs bosons to two gluons are also sizable. For large tan/3 values, i.e. in the 
region of parameter space probed by current and near-future direct Dark-Matter 
search experiments, the total corrections can exceed a factor of three. 



1. Introduction 



There is convincing evidence that dark (non-luminous) matter (DM) contributes between 
20 and 90% of the critical density of the Universe, with values around 30 to 40% currently 
being favored ||T|. Studies of Big Bang nucleosynthesis show that this DM cannot be 
formed by baryons 0, unless they are bound in primordial black holes. Light neutrino DM 
is strongly disfavored by analyses of galaxy formation M . One thus needs physics beyond 
the Standard Model (SM) to explain the dark matter in the Universe. The probably 
best motivated, and certainly most frequently studied, particle physics candidate is the 
lightest neutralino x [@] predicted by supersymmetric extensions of the SM. In the Minimal 
Supersymmetric extension of the Standard Model (MSSM) ||, its stability against decay 
into ordinary particles can be guaranteed by a symmetry called R— parity 0, and for 
wide regions of the supersymmetric parameter space one predicts a relic density in the 
desired range. 

These arguments in favor of neutralino DM have motivated several experimental 
groups around the world to search for ambient relic neutralinos. The strength of the 
expected signal in the two most promising search strategies is directly proportional to 
the neutralino-nucleon scattering cross section c^tv; these are the search for high-energy 
neutrinos originating from the annihilation of neutralinos in the center of the Sun or Earth 
(the so-called "indirect detection") |7|], and the search of the elastic scattering of ambi- 
ent neutralinos off a nucleus in a laboratory detector ("direct search") ||. An accurate 
calculation of u^n for given model parameters is thus essential for the interpretation of 
the results of these searches. 

The matrix element for xN scattering, mediated by Higgs bosons, Z boson and squark 
exchange diagrams, receives both spin-dependent and spin-independent contributions. 
The former are important for neutralino capture in the Sun, but are irrelevant for capture 
in the Earth, and play a subdominant role in most direct search experiments, which 
employ fairly heavy nuclei (Si, Ge, I, Na). The spin-independent contribution in turn 
is usually dominated by Higgs exchange diagrams, where the Higgs bosons couple either 
directly to light (u, d, s) quarks in the nucleon, or couple to two gluons through a loop of 
heavy (c, b, t) quarks or squarks. In some cases squark exchange contributions can also 
be important. 

In this note we point out that these spin-independent contributions receive large 
QCD corrections. The effective Higgs couplings to gluons receive standard QCD cor- 
rections with large positive coefficients, related to the QCD (3— function. In addition, 
gluino-squark loop corrections to the quark masses (or to the Yukawa couplings, for fixed 
physical quark masses) can become very large for down-type quarks, and for the s quark 
in particular. We compute both kinds of corrections, and argue that corrections to other 
xN scattering diagrams should be small. The required formulae are given in Section 2, 
while Section 3 contains some numerical examples. In Section 4 we briefly summarize our 
results and draw some conclusions. 
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2. Formalism 



The spin-independent xN scattering amplitude receives contributions from Higgs bosons 
H and squark |ll[ exchange diagrams. Models with unified gaugino masses typically 
have squark masses much larger than the neutralino mass, rrtg £ 5m^, in which case the 
dominant contribution usually comes from Higgs boson exchange. Note that only scalar 
Higgs couplings to neutralinos contribute in the non-relativistic limit. In the absence of 
significant CP-violation in the Higgs sector, one therefore only has to include contributions 
of the two neutral CP-even Higgs bosons. The contribution of the heavier Higgs boson 
often dominates, since its couplings to down-type quarks are enhanced if the ratio of 
vacuum expectation values (vev's) is large, tan/5 ^> 1. 

CP-even Higgs bosons can couple to nucleons in two different ways: through an ef- 
fective coupling to two gluons mediated by loops of heavy quarks [12| or squarks |Tl[ 



or 



directly to light quarks in the nucleon. The leading contribution to the Hfgg couplings 
comes from heavy quark triangle diagrams. To leading order, they can be described by 
the effective Lagrangian 
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where H® is one of the neutral CP-even Higgs bosons h and H and F^ ua is the gluon field 
strength tensor, a being a color index. Note that the dimensionless coefficients qq are 
independent of mg; the factor mg in the H®QQ coupling is canceled by a factor 1/rriQ 

Explicit expressions for these coefficients can e.g. be found 
The effective Lagrangian eq. (ffl) gives rise to HfNN couplings, N being a 



from the loop integral |12 
in Ref. |TIJ 



nucleon, through hadronic matrix elements [|TT|, |T!| 
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The first point we wish to make is that the effective interaction ([!]) is subject to large 
corrections from ordinary QCD. These can most easily be computed using low-energy 
theorems 
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The effective Higgs-gluon interaction can then be rewritten as 
1 Pq(®s) r nuvaSr Ci Q 
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Here [3q is the contribution of heavy quark Q to the QCD f3— function and 7q is the 
anomalous dimension of rriQ. Eq.(|]) only describes interactions where all three external 
legs couple to the heavy quark line, but is valid to all orders in perturbation theory. 
Currently results up to order are known 
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The scale of the overall factor a s need not be specified, since the matrix element (0) is 
scale-independent. However, the size of the higher order corrections in the square brackets 
are determined by a s taken at the scale rriQ, which is the only energy scale in the problem 
(note that we are interested in configurations with essentially vanishing momentum flow 
through the diagram). These corrections are therefore larger for Q = c than for Q = t. 
Nf counts the number of active flavors at scale tuq; we take Nf = 3, 4, 5 for Q = c,b, t. 

The general result eq. can also be used for squark loop contributions to the Higgs- 
gluon coupling. Here results up to 0(a 2 s ) are known, giving a correction factor -|> "^"^ 



in the effective Lagrangian, relative to the leading order result listed in Ref. [IT|. Note that 
the correction factor is even larger than for quark loops; however, the overall contributions 
of squark loops to the effective H®gg couplings at vanishing external momenta are always 
much smaller |Tl| than the quark loop contributions. 

Potentially even larger corrections come from gluino-squark loop contributions to 
the H®qq couplings, which are closely related to g — q loop corrections to quark masses 
[T5|j . These loop contributions induce a coupling of quarks to the "wrong" Higgs current 
eigenstate. In particular, the masses of down-type quarks will now receive contributions 
from the vev of the Higgs field with positive hypercharge. Parametrically, these corrections 
are therefore of the order 5md/md ~ —tan/3. These corrections can thus become 0(1) 
for tan/3 ~ 50, which is within the allowed range. 

The large size of these corrections raises concerns regarding the reliability of pertur- 



bation theory. However, it has recently been shown |T6|, [17J that the one-loop corrections 
can be written in such a way that all higher-order corrections of the form (a s /ir tan/3) n 
are re-summed. Note that the corresponding corrections to the couplings of up-type 
quarks are much smaller for tan/5 > 1; given the uncertainties of the hadronic matrix 
elements appearing in the xiV scattering amplitude, corrections of order 10% or less can 
safely be ignored. We are thus only interested in the couplings of down-type quarks D 
to CP-even Higgs bosons h, H, which can be written as [17] 



4(M S ) = f-r4 M (M s ) cos(a - 0) + m^ SSM (M s ) sin(a - 0) tan/?! (6) 
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s°(M s ) = — 9 — \mjf(M s ) sin(a - 0) + m^ SSM (M s ) cos(a - 0) tan/31 (7) 



2M, 



where g is the SU(2) gauge coupling, and a is the Higgs mixing angle in the notation of 
Ref. . In eqs. (0) m| M (M s ) stands for the running quark masses, where only standard 
QCD corrections are included; their numerical values are known (with some uncertainties) 
from experiment. m^ SSM (Ms) is just the product of the D Yukawa coupling with the 
vev of the Higgs field with negative hypercharge; it differs from nr^ 1 by sparticle loop 
corrections: 



m^ SSM (M s ) = m^ M (Ms) - ^^m 5 sin(2^) [B (ml,m~ 9 ,m 3l ) - B^m^m^ 



(8) 
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Here 9 D is the Dl — Dr mixing angle, m Dl and m D are the smaller and larger eigenvalue 
of the D squark mass matrix, rrig is the gluino mass, and Ms is a scale of the order of the 

squark or gluino mass; we will take the value Mg = lmjm Di m D2 J . The first argument 
of the Passarino-Veltman two-point function Bq can, to excellent approximation, be set 
to zero. In the case of the bottom quark (which contributes to the H®gg couplings) we 
add an additional correction from chargino-stop loops, which comes from the large top 
Yukawa coupling: 
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Note that for tan/? 3> 1, sm.(26 D ) ~ 2m D jj,tan(3/(m 2 ^ — rri^), and 1/ sin(2/9) ~ tan/5/2; 
the corrections eqs. (§) and (|) are therefore enhanced for large tan/3, as advertised. Of 
course, if these corrections are set to zero, eqs. (0) reduce to the standard tree-level H®qq 
couplings 



It is well known that Yukawa couplings of quarks receive large QCD corrections, i.e. 
they show a strong scale dependence. Here we need these couplings at a relatively low 
scale, of the order of the quark mass itself, or in case of light quarks, at a scale not far from 
1 GeV. In contrast, eqs. (|7|) should be interpreted as boundary conditions of the couplings 
of the effective theory that emerges when squarks and gluinos are integrated out; these 
boundary conditions are thus valid at a scale Ms. Below this scale these couplings run 
according to standard renormalization group equations (RGE) |Tj| : 
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where Nf is again the number of active flavors, and Qi and Q2 are two energy scales. 
Since this number changes whenever we cross a heavy quark threshold, the total effect of 
the running has to be computed by applying eq. (|T0"1) repeatedly; e.g., 



h (M s 



'a a (M s )~ 


12/21 


a s {m t ) 


12/23 


a s (m b ) 


a s (m t ) _ 




_a s (m b )_ 




_a s (m c )_ 



-1 12/25 



(11) 



We have seen in eq. (^) that the size of the g — q loop corrections to mo is proportional 
to the amount of Dl — Dr mixing, which in turn is proportional to in o- This raises the 
question which m D to use when computing sin(26 ) £ ) ). Since the corrections eqs. (^) and (^) 
are to be applied at the high scale Ms it is fairly clear that one should use a running mass 
mo at this high scale. At the one-loop order we cannot distinguish between the choices 
here. However, in Ref. [|1^] it has been shown that loop corrections 



m D M and m D 



MSSM 



to Higgs boson decay widths into quarks and squarks will be small only if the tree- 
level squark mass matrices are written in terms of m^ SSM , i.e. if m|J SSM is used in 
the calculation of sin(2#£,). This requires an iteration, which however usually converges 
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quickly. This iteration is equivalent to the re-summation of higher orders advocated in 
Ref. | I7[l . We will show below that the difference between using mf) M and m^ SSM in the 



calculation of sin(26^) can indeed be numerically significant. 

The squark-gluino loop corrections to the mass of d— type quarks also affect the lead- 
ing 0{rriz 2 ) spin-independent contributions from d— type squark exchange. These con- 



tributions are proportional to the quark masses [TO], |Tf[, either through the interference 
of gauge and Yukawa contributions to the neutralino-quark-squark couplings (which re- 
quires gaugino-higgsino mixing in the neutralino sector), or through q^ — q~R mixing. 
These corrections can again most easily be understood in terms of an effective f q qqxx 
interaction, where the coefficient f q is determined by matching to the full theory at a 
scale Q ~ m q [TT|] . Since both the (Jl — qn mixing angle and the quark Yukawa coupling 
are determined by m^ 18 ^ rather than by m^ M , the gluino-squark loop correction to the 
effective f q qqxx interaction can simply be obtained by multiplying the usual tree-level 
result |T(], [TIJ] with m^ lsSM (Ms)/m q iM (M s ). Of course, we need these couplings at a low, 
hadronic scale. However, the scale dependence of these couplings is identical to that of 
the Hqq couplings, i.e. is described by eq. ([[(]). This amounts to a purely multiplicative 
renormalization, which leaves the relative size of the g — q loop corrections unchanged. 
Altogether we thus have 

-vnMSSMf T\/f \ 

AD), _ ad), rn D [Ms) 

JD I improved ~ JD I tree " SM/i/ \ ' 

m D { M S) 

where the superscript (D) signifies that this expression only applies to the D squark 
exchange contribution to the effective DDxx interaction^. 

The remaining potentially important contributions to the xiV scattering amplitude 
are due to gauge interactions, at least in the more appealing scenario where the LSP is 
Bino- or photino-likeQ These include spin-dependent contributions due to Z and squark 



exchange, as well as spin-independent 0(rrig ) contributions fil l. Since to one-loop 



order electroweak gauge couplings are not renormalized by strong interactions, we do not 
expect significant QCD corrections to these contributions. However, as already stated in 
the Introduction, these two contributions are often subdominant. 



3. Results 

We are now ready to present some numerical results. In order to illustrate the importance 
of the corrections described in the previous Section, we show examples for the ratio R of 
the neutralino scattering rate on 76 Ge with and without these corrections. If the small 

1 In fact, in the "decoupling limit", where cos(a — j3) — ► 0, the total g — q loop correction to the H 
exchange contribution takes the same simple form as in eq. (|l2l), while the h exchange contribution is 
not affected by these corrections. 

2 We consider this to be more appealing since one then can obtain an interesting relic density with 
sparticle masses in the "natural" domain of a few hundred GeV. Higgsino-like LSPs need masses around 
a TeV or more to contribute significantly to the mass density of the Universe. 
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difference between the x n an d XV scattering amplitudes is neglected, R is simply the ratio 
of the corrected and uncorrected %iV scattering cross sections. 

In Fig. (la,b) we show results for a scenario where all soft SUSY-breaking contribu- 
tions to sfermion masses at the weak scale have a common value mo = 600 GeV. Similarly, 
we took a common value A = 1.2 TeV for all trilinear soft SUSY-breaking parameters at 
the weak scale. We chose a relatively small mass of the CP-odd Higgs bosons, = 240 
GeV ~ 1.6m^, in order to obtain LSP detection rates of interest for present experiments, 
at least in the region of large tan/3. We assumed the usual unification conditions for 
gaugino masses, with an SU (2) gaugino mass M 2 = 300 GeV at the weak scale; note that 
now also the gluino mass becomes relevant, see eq. (H). We present results for = M 2 
(Fig. la) and = 2M 2 (lb), as a function of tan/3. In both cases we show results for 
both positive (dashed) and negative (solid) fi. Finally, the upper (lower) curve with a 
given pattern has been obtained using (m^ M ) when calculating sin(2^). 



M 2 = 300 GeV, m =2.5m A =A /2 = 600 GeV 




Q g u i i I i i i I i i i I i Ll I Q Li i i I i i i I i i i I i i i 

20 40 60 80 J 20 40 60 80 
tan/? tan/3 

Figure 1: The ratio R of corrected to uncorrected \N scattering cross section in a model 
with universal soft breaking parameters at the weak scale. The upper (lower) curve 
of a given pattern uses the quark mass with (without) sparticle loop corrections when 
computing the squark L-R mixing angle. 
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For small values of tan/3 the total corrections are dominated by the standard QCD 
corrections (|5[) to the H®gg couplings. They increase this contribution to the total X-/V 
scattering amplitude by ~ 25%, which leads to an increase of the counting rate by ~ 15%. 

As anticipated the g — q loop corrections can be much larger, if tan/5 ^> 1. Our 



sign convention for \i coincides with that of Refs. Jy], |18[. For tan/3 ^> 1 the L — R 
mixing of d— type squarks is essentially proportional to /i; we then find a negative g — q 
loop correction to mp, and hence a positive correction to the Yukawa coupling for fixed 
physical quark mass, if /x < 0. Since it is the Yukawa coupling, rather than the quark 
mass, which appears in most spin-independent contributions to the scattering amplitude, 
/i < therefore also leads to positive corrections to the scattering cross sectionP]. Not 
surprisingly, these corrections are considerably larger for larger |/i|, as can be seen by 
comparing Fig. la with Fig. lb. 

Note that we have (somewhat arbitrarily) terminated the curves when the bottom 
Yukawa coupling hf,(Ms) > 1.2, which roughly corresponds to its "fixed point" value. An 
even larger hb{Ms) would lead to a Landau pole at higher energies, i.e. the Yukawa sector 
of the theory would become non-perturbative at energies only slightly above Ms. Since 
for a fixed value of tan/3 the g — q corrections increase the Yukawa coupling if /j < 0, 
for this sign of /i the upper bound on tan/? becomes significantly stronger; conversely, 
for fi > larger values of tan/3 become accessible once the corrections are included. 
Since the predicted counting rate increases with increasing tan/3, the maximal possible 
\N scattering cross section for fixed values of the dimensionful input parameters, which 
occurs for the biggest allowed value of tan/3, would thus not be affected by the corrections 
if the relative size of these corrections were the same for d, s and b quarks. 

Even though we have assumed flavor-independent soft breaking parameters, this con- 
dition is not satisfied. First, while L-R mixing does not change the masses of the d and 
s squarks very muchF], it does affect the b masses significantly. The corrections (||) are 
therefore larger, as fraction of the tree-level mass, for D = b than for D = d,s. Moreover, 
the chargino-stop loop correction (||) only contributes to the b— quark mass. Note that it 
has the opposite sign as the g — q loop correction. For the parameters of Fig. 1 this is the 
dominant effect. In fact, in Fig. la this top Yukawa correction is nearly as large as the 
SUSY-QCD correction. The reason is that ti — tn mixing, which is relevant in eq. @, is 
proportional to A t (= Aq in our case) which is four times bigger than \fi\ in this example. 
As a result, for /i < the loop corrections increase the predicted counting rate at the 
largest allowed value of tan/3 by a factor of 1.9 (1.4) in Fig. la (lb). On the other hand, 
the increase of the maximal allowed value of tan/3 for /i > implies that for this sign of /i 
the maximal rate only decreases by a factor 0.79 (0.74) in Fig. la (lb); of course, one also 
has to keep in mind that the corrections (|5]) to the H®gg couplings are always positive. 

3 The cross-over of the two sets of curves at small tan/3 in Fig. la occurs because at small tan/3 and 
/i > the most important contribution comes from the h boson exchange, whereas for /i < the H boson 
exchange is always dominant. 

4 The difference between the squark masses is in these cases mostly determined by the different D— term 
contributions to the masses of SU(2) singlet and doublet squarks. 



7 



Note that using m^ SSM rather than m 8M in the calculation of sin(2^) always increases 
the prediction of the loop-corrected xN scattering rate. Recall that a negative correction 
to the quark mass implies an increase of the scattering cross section. Since in this case 
m MSSM > rriq M , using 77i^ ISSM increases L-R mixing and hence gives a larger, positive 
correction to the predicted rate. Conversely, if the corrections to the quark masses for 
fixed Yukawa couplings are positive, the correction to the xN scattering amplitude is 
negative. Since now 77i^ SSM < m 8M , the use of m^ 1381 ^ reduces L-R mixing, and hence 
lowers the absolute size of the negative corrections to the scattering rate. 

The corrections to hb for fixed physical m& can have even larger effects if one as- 
sumes universality of soft breaking parameters at some very high scale, as in the popular 
mSUGRA models with radiative SU(2) x U(l)y symmetry breaking |J. Here one as- 
sumes a common squared soft breaking mass Wq not only for sfermions, but also for 
Higgs bosons; however, this universality only holds at the scale of Grand Unification, 
M x ^ 2 • 10 16 GeV. The squared soft breaking mass of one of the Higgs bosons is then 
driven to negative values by the top Yukawa contribution to the relevant RGE. One can 
show [^(J that the physical mass of the CP-odd neutral Higgs boson is then given byf] 



m\ = ^f + 0(ht), (13) 
sin lj 

where the corrections from the bottom Yukawa coupling hb are negative. In fact, the 
upper bound on tan/3 often (but not always) comes from the requirement that itla should 
be sufficiently large. Note that the relevant quantity here is again the bottom Yukawa 
coupling, the size of which depends on the size and sign of the g — q and t — x^ 1 loop 
corrections. The value of uia is in turn strongly correlated with the mass tuh of the 
heavier neutral CP-even Higgs boson, the exchange of which usually dominates the spin- 
independent xN scattering amplitude for tan/3 3> 1. In this kind of model one thus expects 
a much stronger dependence of the predicted xiV scattering rate on these corrections than 
in models where all soft breaking parameters are fixed at the weak scale, independently 
of tan/3. 

This is illustrated in Fig. 2, where we took a weak-scale £77(2) gaugino mass of 200 
GeV, and GUT-scale boundary conditions mo = Aq = 250 GeV. Note that in this model 
the absolute size of fi is a derived quantity, but its sign can still be chosen freely; dashed 
(solid) lines are again for positive (negative) \i. The two inner curves have been obtained 
by using m? M when calculating L-R squark mixing, as well as hb as it appears in the 
RGE. In this case the correction "only" amounts to a factor of 2 at most, similar to the 
case shown in Fig. la. 

The outer curves in Fig. 2 have been obtained by using both in the calculation 

of sm(29q) and in the calculation of the bottom quark Yukawa coupling in the RGE. 
Mostly due to the change of tua discussed above, we now find correction factors as large 
as 100 for fixed values of tan/3. Note that in the given example for fi < the upper bound 

5 We ignore small weak-scale threshold corrections. 
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mSUGRA: M 2 =200 GeV, m =A = 250 GeV 



0.1 



0.01 




tan/3 



Figure 2: The ratio R of corrected to uncorrected xN scattering cross section in an 
mSUGRA model. The outer (inner) pair of curves uses the quark mass with (without) 
sparticle loop corrections when computing the squark L — R mixing angle as well as the 
6-quark Yukawa contribution to the RGE. 



on tan/3 is determined by the lower bound on m^; we have required rriA + mf l > 180 (150) 
GeV if cos 2 (a — (3) ~ 1 (> 0.2), to cope with the experimental constraints. As a result 
the predicted counting rate at the maximal allowed value of tan/3 changes very little when 
the loop corrections are included.^ On the other hand, for /i > the upper bound on 
tan/3 after the inclusion of sparticle loop corrections comes from the requirement that the 
lighter f mass eigenstate should not be lighter than the lightest neutralino. In this case 
the corrections therefore reduce the maximal possible b (and s and d) Yukawa couplings, 
leading to a reduction of the predicted \N scattering rate at the maximal allowed value 
of tan/3 by a factor close to 100. 

6 Note that we always include the g — q loop corrections to the mass matrix of neutral Higgs bosons 
PH . The dominant contribution here usually comes from i — g loop corrections to rrit , but for large tan/3 
and small ttia values, b — g loop corrections to mb can also have a significant effect on the Higgs mixing 
angle a. 
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4. Summary and Conclusions 



In this paper we discussed two kinds of (potentially) large QCD corrections to the xN 
scattering amplitude. Standard QCD corrections increase the effective coupling of CP- 
even Higgs bosons to two gluons, which proceeds through a loop of heavy quarks, by ~ 
25%. Supersymmetric QCD corrections to the quark masses, or to their Yukawa couplings 
for fixed physical quark masses, can have even larger effect if tan/? 3> 1: we found that 
an enhancement or suppression of the predicted xN scattering cross section by a factor 
of 2 is quite easily possible for a fixed set of soft breaking parameters and fixed tan/3. 
Since the correction can also change the maximal allowed value of tan/?, the change of 
the maximal possible counting rate, which occurs at the largest allowed value of tan/3, is 
usually smaller, but can still be significant. 

Much bigger effects from g — q loop corrections are possible if the size of the b Yukawa 
coupling affects the particle spectrum at the weak scale, as in mSUGRA models. The 
dominant effect here comes from the anti-correlation between the b Yukawa coupling and 
the masses of the heavier Higgs bosons. Since for large tan/3 the scattering cross section 
roughly scales like the inverse fourth power of the mass of the heavier neutral CP-even 
Higgs boson, the sparticle loop-induced change of the b Yukawa coupling can change the 
predicted neutralino detection rate by up to two orders of magnitude, if tan/3 is large .[] 
Note that tan/3 3> 1 is required in this model to obtain scattering rates that are sufficiently 
large to be probed in present or near-future experiments. 

One consequence of these sparticle loop corrections is a strong dependence of the 
predicted xN scattering rate on the sign of [i even for tan/3 3> 1. It is well known that 
for small and moderate tan/3, positive values of \i lead to more gaugino-higgsino mixing 
in the neutralino eigenstate, and hence to larger couplings of the neutralino to Higgs 
bosons, than do negative values of \i. This tree-level dependence of the scattering rate 
on the sign of \i disappears at large tan/3. However, there the g — q loop corrections 
become important. Note that now fi < gives larger counting rate (for fixed than 
/i > 0. These corrections can thus qualitatively change the dependence of the predicted 
relic neutralino detection rate on the parameters of the supersymmetric model. 
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